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I.  INTRODUCTION 


Firing  tests1  of  the  PXR  projectiles  with  modified  base  geometries  have 
recently  been  performed.  Analysis  of  the  test  results  revealed  small,  but 
experimentally  significant,  differences  in  the  observed  effective  drag. 
Therefore,  it  was  desired  to  numerically  compute  the  flow  field  for  these 
configurations  and  determine  the  effect  of  the  modified  base  geometries  on  the 
base  region  flow.  The  total  drag  for  projectiles  can  be  divided  into  three 
components:  (i)  pressure  drag  (excluding  the  base  region),  (ii)  viscous  or 
skin  friction  drag,  and  (ill)  base  drag.  At  transonic  speeds  the  base  drag 
constitutes  a  major  portion  of  the  total  drag.  Thus,  the  determination  of 
base  pressure  for  projectile  configurations  is  essential  in  predicting  the 
total  drag  coefficient. 

A  series  of  computations  have  recently  been  performed  for  a  PXR  projec¬ 
tile  configuration  in  which  the  full  range  of  aerodynamic  coefficients  have 
been  determined.2  The  actual  PXR  configuration  is  a  rocket  assisted  projec¬ 
tile  and  contains  a  base  cavity  as  part  of  the  base  configuration.  The  pre¬ 
vious  computations2  however,  did  not  include  the  cavity  and  modeled  the  base 
as  a  solid  wall.  In  pursuing  the  capability  to  compute  the  flow  field  of  the 
actual  projectile,  a  new  capability  has  been  developed  which  allows  the  base 
to  include  an  annular  cavity.  Other  base  conf i gurati ons  which  contain  cylin¬ 
drical  cavities  have  not  been  considered  in  the  present  analysis.  However, 
with  minor  modifications  the  developed  code  could  handle  such  cavities.  This 
report  discusses  the  development  of  this  capability  and  the  preliminary 
results  which  have  been  obtained. 


II.  COMPUTATIONAL  TECHNIQUE 

The  Azimuthal  Invariant  (or  Generalized  Axi symmetri c)  thin-layer  Navier- 
Stokes  equations  for  general  spatial  coordinates  £,  n,  C  can  be  written  as3 


3  q  +  51  +  3G  +  H  =  Re18S 


(1) 


1 .  Private  Communications  with  Mr .  Richard  Eitemiller  of  the  Firing  Tables 
Branch ,  Launch  and  Flight  Division ,  BRL,  June  1983 . 

2.  C.  J •  Nietubicz,  R .  LaFarge,  J .  Sahuy  and  I) .  C •  My  tin,  "Aerodynamic  Coef¬ 
ficient  Predictions  for  a  Projectile  Configuration  at  Transonic  Speeds 
AIAA  Paper  No*  84-0326,  January  1984. 

3.  C.  J .  Nietubicz,  T.  H.  Pulliam ,  and  J .  L.  Steger,  "Numerical  Solution  of 
the  Azimuthal-Invariant  Thin-Layer  Navier-Stok.es  Equations JJ.S.  Army 
Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground ,  Maryland ,  ARBRL - 
TR-02227,  March  1980.  (AD  A085716)  (Also  see  AIAA  Paper  No.  79-0010, 
January  1979. 
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where 


?  =  ?(x,y,z,t)  is  the  longitudinal  coordinate 
o  =  n(y,z,t)  is  the  circumferential  coordinate 


C  =  c(x,y,z,t)  is  the  normal  coordinate 
x  =  t  is  the  time 
and 
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represent  the  contravariant  velocity  components.  Equation  (1)  is  solved  in  a 
time  asymptotic  fashion  with  interest  only  in  the  steady-state  solution.  The 
numerical  algorithm  used  is  a  fully  implicit,  approximately  factored  finite 
difference  scheme.  The  algorithm  is  first  order  accurate  in  time  and  fourth 
order  in  space.  A  two-layer  algebraic  eddy  viscosity  model1*  is  included  for 
the  computation  of  turbulent  flows.  Details  of  the  assumptions  and  the 
algorithm  are  included  in  References  5-7. 


III.  METHOD  OF  SOLUTION 
A.  Base  Region  Flow  without  Cavity 

The  procedure  used  to  compute  the  base  flow  for  shell  without  the  base 
cavity  has  been  described  in  Reference  8;  however,  limited  details  are 
repeated  here  for  clarity.  The  base  flow  code  computes  the  full  flow  field 
(including  the  base  region)  of  a  projectile.  Figure  1  shows  a  schematic 
illustration  of  the  flow  field  segmentation  and  shows  the  transformation  of 
the  physical  domain  into  the  computational  domain.  This  flow  field  segmenta¬ 
tion  procedure  is  equivalent  to  using  multiple  adjoining  grids.  An  important 
advantage  of  this  procedure  lies  in  the  preservation  of  the  sharp  corner  at 
the  base  and  allows  easy  blending  of  the  computational  meshes  between  the 
regions  ABCD  and  AEFG.  No  approximation  of  the  actual  sharp  corner  at  the 
base  is  made.  Thus,  realistic  representation  of  the  base  is  inherent  in  the 
current  procedure. 

The  cross  hatched  region  in  the  physical  domain  represents  the  model. 
The  line  BC  is  the  base  and  the  region  ABCD  is  the  base  region  or  the  wake. 
The  line  AB  is  a  computational  cut  through  the  physical  wake  region  which  acts 
as  a  repetitive  boundary  in  the  computational  domain.  Implicit  integration  is 
carried  out  in  both  5  and  c  directions  (see  Figure  1).  Note  the  presence  of 


4.  B.  S.  Baldwin,  and  H.  Lomax,  " Thin-Layer  Approximation  and  Algebraic  Model 
for  Separated  Turbulent  Flews,"  AIAA  Paper  No.  78-257,  1978. 

5.  R.  Beam,  and  R.  F.  Warming,  "An  Implicit  Factored  Scheme  for  the  Compres¬ 
sible  Navier-Stdk.es  Equations AIAA  Journal ,  Vol.  16,  No.  4,  April  1978, 

pp.  393-402. 

6.  J.  L.  Steger,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoils,"  AIAA  Journal ,  Vol.  16, 
No.  4,  July  1978,  pp.  679-686. 

7 .  T.  H.  Pulliam,  and  J.  L.  Steger,  "On  Implicit  Finite-Difference  Simula¬ 
tions  of  Three-Dimensional  Flow,"  AIAA  Jourmal ,  Vol.  18,  No.  2,  February 
1980,  pp.  159-167. 

8.  J.  Sahu,  C.  J.  Nietubicz,  and  J.  L.  Steger,  "Numerical  Computation  of  Base 
Flow  for  a  Projectile  at  Transonic  Speeds,"  U.S.  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02495 ,  June  1983. 
(AD  A130293)  (Also  see  AIAA  Paper  No.  82-1358,  August  1982). 
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the  lines  BC  (base)  and  EF  (nose  axis)  in  the  computational  domain.  They  both 
act  as  boundaries  in  the  computational  domain  and  special  care  must  be  taken 
in  inverting  the  block  tridiagonal  matrix  in  the  5  direction.  The  details  of 
these  can  be  found  in  References  8  and  9  and  are  not  included  here. 

The  no  slip  boundary  condition  for  viscous  flow  is  enforced  by  setting 


U  =  V  =  W  =  0 


(2) 


on  the  body  surface.  At  the  base  boundary,  an  inviscid  boundary  condition  has 
been  used.  Additionally,  at  the  corner  of  the  base,  the  boundary  conditions 
are  double-valued  and  depend  on  the  direction  from  which  the  corner  is 
approached.  Approaching  the  corner  in  the  streamwise  direction  (EB),  the  no¬ 
slip  boundary  condition  is  used  while  approaching  in  the  radial  direction 
(along  the  base,  CB) ,  the  inviscid  boundary  condition  is  used.  Along  the 
computational  cut  (AB) ,  the  flow  variables  above  and  below  the  cut  are  simply 
averaged  to  determine  the  boundary  conditions  on  the  cut.  On  the  centerline 
of  the  wake  region,  a  symmetry  condition  is  imposed  and  free  stream  conditions 
are  used  on  the  outer  boundary. 

B.  Base  Region  Flow  with  Cavity 

The  procedure  used  to  compute  the  base  flow  including  the  effect  of  the 
cavity  is  now  described.  The  flow  field  segmentation  procedure  described 
above  takes  a  new  form  shown  schematically  in  Figure  2.  There  are  two  flow 
field  regions  AEFG  and  ABPRSQCD.  The  latter  one  contains  the  base  cavity 
PRSQ.  The  two  cross  hatched  regions  in  the  computational  domain  are  the  aft 
end  of  the  projectile  and  computations  are  not  made  in  these  regions.  The 
presence  of  such  rectangular  regions  inside  the  computational  domain  requires 
special  care  in  carrying  out  the  integrations  both  in  the  5  and  the  c  direc¬ 
tions.  Modeling  the  base  cavity  in  this  manner  requires  specification  of 
boundary  conditions  on  the  new  boundaries  PR,  RS  and  SQ  and  special  monitoring 
of  the  5  and  ?  indices. 

The  boundary  conditions  on  the  body  surface,  the  cut  and  the  downstream 
boundary  remain  the  same  as  described  earlier.  Along  the  base  (lines  BP,  QC 
and  RS)  inviscid  boundary  conditions  are  used.  Inviscid  boundary  conditions 
are  also  used  for  lines  PR  and  SQ.  The  turbulence  model  has  not  been  modified 
for  this  region  and  the  same  length  scales  used  for  the  near  wake  are  used  in 
the  cavity. 


9. 


J.  Sahu,  C.  J.  Nietubicz,  and 
Projectile  Base  Flow  with,  and 
Research  Laboratory,  Aberdeen 
November  1983.  (AD  A135783) 
1983)  . 


J.  L.  Steger,  "Navier-Stokes  Computations  of 
without  Base  Injection,"  U.S.  Army  Ballistic 
Proving  Ground,  Maryland,  ARBRL-TR-02532, 
(Also  see  AIAA  Paper  No.  83-0224,  January 
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IV.  MODEL  AND  GRID 


The  basic  model  geometry  used  in  this  study,  excluding  the  details  of  the 
base  is  shown  in  Figure  3.  It  is  a  secant-ogi ve-cyl inder-boattai 1  configura¬ 
tion.  The  base  region  configurations  of  interest  are  shown  in  Figure  4. 
These  configurations  all  contain  annular  cavities.  Free  flight  tests  have 
previously  been  conducted  on  a  SAMOS  re-entry  configuration10  which  contained 
a  cylindrical  cavity.  Cylindrical  cavities  are  not  included  in  the  present 
analysis.  In  this  report,  the  annular  base  cavity  shown  in  Figure  4(a)  has 
been  modeled. 

The  computational  grid  used  for  the  numerical  calculations  was  obtained 
from  a  grid  generator  described  in  Reference  11.  This  program  allows  arbi¬ 
trary  grid  point  clustering,  thus  enabling  grid  points  for  the  projectiles  to 
be  clustered  in  the  vicinity  of  the  body  surface.  The  full  grid  is  shown  in 
Figure  5  while  Figure  6  shows  an  expanded  view  of  the  grid  near  the  model. 
The  grid  consists  of  114  points  in  the  longitudinal  direction  and  40  points  in 
the  normal  direction.  For  the  base  region,  30  grid  points  were  used  in  the 
streamwise  direction.  The  grid  nodes  in  the  normal  direction  are  clustered 
near  the  surface  to  capture  viscous  effects.  These  points  were  exponentially 
stretched  away  from  the  surface  with  a  minimum  spacing  of  .00002D  at  the  wall. 
This  locates  at  least  two  points  in  the  laminar  sublayer.  As  can  be  seen  in 
Figure  6,  the  grid  points  in  the  longitudinal  direction  are  clustered  near  the 
cylinder  and  boattail  junctions  and  also  at  the  base  where  large  gradients  in 
flow  variables  are  expected  to  occur. 

An  expanded  view  of  the  grid  in  the  base  region  for  the  case  of  a  solid 
base  is  shown  in  Figure  7.  Figure  8  shows  the  grid  in  the  base  region  includ¬ 
ing  the  cavity.  For  computations  of  base  flow  with  base  cavity  it  is  neces¬ 
sary  to  locate  an  adequate  number  of  grid  points  in  the  cavity  region  and  the 
exponential  stretching  used  previously  is  not  suitable.  Therefore,  a  1-D 
elliptic  solver  was  used  to  obtain  the  grid  spacings  in  the  normal  direction  l 
for  the  base  region  with  the  cavity.  In  both  cases  with  and  without  the 
cavity,  the  minimum  spacing  on  both  sides  of  the  cut  (AB)  was  the  same 
(.00002D).  Thus,  a  smooth  variation  of  the  grid  is  maintained  across  the  cut. 


V.  RESULTS 

Computations  for  the  PXR  projectile,  both  with  and  without  the  cavity, 
have  been  made  at  M  =  0.9  and  a  =  0.  Both  qualitative  and  quantitative 
results  are  now  presented  for  these  base  configurations.  Figure  9  shows  the 


10.  E.  D.  Boyer,  "Free  Flight  Tests  of  the  SAMOS  Re-entry  Configuration ," 
U.S.  Army  Ballistic  Research  Laboratory,  Maryland,  ARBRL-MR-1321 , 
February  1961.  (AD  374655L) 

11.  J.  L.  Steger,  C.  J.  Nietubicz,  and  K.  R.  Heavey,  " A  General  Curvilinear 
Grid  Generation  Program  for  Projectile  Configurations U.S.  Army 
Ballistic  Research  Laboratory,  Maryland,  ARBRL-MR-03142,  October  1981. 
(AD  A1 07334) 
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velocity  vectors  in  the  base  region  for  the  case  with  no  cavity.  The  recircu- 
latory  flow  in  the  near  wake  is  clearly  evident.  The  reattachment  point  is 
less  than  a  caliber  downstream  of  the  base.  The  effect  of  the  cavity  on  the 
base  region  flow  field  is  shown  in  Figure  10.  The  base  region  flow  field  is 
now  characteri zed  by  two  recirculation  regions  of  opposite  sign.  A  small 
bubble  is  formed  inside  the  cavity  which  has  displaced  the  shear  layer.  The 
reci rcul ati on  bubble  in  the  near  wake  is  now  larger  and  the  reattachment  point 
has  been  moved  further  downstream  to  a  position  1.5  calibers  from  the  base. 

A  more  critical  look  at  the  computational  results  is  presented  in  Figures 
11  and  12.  Recall  that  the  computational  technique  used  in  these  calculations 
is  a  time  dependent,  implicit  algorithm  and  the  unsteady  set  of  equations  are 
marched  in  time  to  obtain  the  steady  state  results.  An  implicit  algorithm,  in 
general,  allows  large  time  steps  in  order  to  reach  the  desired  steady  state 
solution.  However,  large  time  steps  can  often  give  rise  to  numerical  oscilla¬ 
tions  as  evidenced  by  results  shown  in  Figures  11  and  12.  Figure  11  and  12 
show  the  variation  of  base  drag  and  total  drag  respectively  with  time.  With  a 
non-dimensional  time,  At  =  .01,  the  solutions  indicate  an  unsteady  flow  field. 
The  time  step  was  reduced  to  .005  and  the  amplitude  of  the  oscillations  drop¬ 
ped.  The  amplitude  of  the  oscillations  was  further  reduced  for  a  time  step  of 
.0025.  This  indicates  that  the  flow  field  in  the  base  region  with  the  cavity 
may  not  be  unsteady.  However,  it  requires  a  different  approach  for  advancing 
in  time  to  get  the  desired  steady  state  results.  This  type  of  oscillation  was 
not  found  in  the  computations  of  the  base  flow  without  the  cavity. 

Such  oscillations  are  believed  to  be  numerical  rather  than  physical  and 
are  dependent  on  the  flow  problem.  Similar  oscillations  were  also  observed  in 
a  time  dependent  numerical  computation12  of  supersonic  flow  over  a  ramp. 
Additionally,  these  oscillations  were  not  observed  when  the  eddy  viscosity 
turbulence  model  was  frozen.  12  This  again  suggests  that  the  oscillations  in 
the  present  calculations  are  not  physical  and  are  probably  due  to  the  inter¬ 
action  of  the  turbulence  model  with  the  flow  field. 

The  total  drag  as  shown  in  Figure  12  is  converging  to  a  lower  value  with 
the  cavity  as  compared  to  the  correspondi ng  value  of  total  drag  obtained  for 
the  solid  base.  Other  drag  components  such  as  pressure  drag  and  skin  friction 
drag  remained  practically  the  same  both  with  or  without  the  cavity.  The 
influence  of  the  cavity  on  the  flow  field  upstream  of  the  base  corner  is  neg¬ 
ligible.  Thus,  no  significant  change  in  pressure  drag  or  skin  friction  drag 
was  found.  The  cavity,  however,  has  a  large  influence  on  the  flow  field  in 
the  base  region,  especially  in  the  near  wake.  Thus,  the  base  drag  is  strongly 
affected.  Any  change  in  the  base  drag  is,  of  course,  reflected  in  the  total 
aerodynamic  drag. 

As  can  be  seen  from  Figure  12,  a  reduction  in  the  total  drag  due  to  the 
cavity  of  about  25%  is  predicted.  The  effect  of  a  base  region  cavity  on  the 


12.  D.  Degani,  and  J.  L.  Stegery  " Comparison  Between  Navier^Stokes  and  Thin- 
Layer*  Computations  for  Separated  Supersonic  FZou)9"  AT  A  A  J  ournaly  Vot .  21y 
No .  lly  November  1983 y  pp.  1604-1606 . 
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base  pressure  has  been  studied  by  other  researchers. 13  The  trend  in  the 
reduction  of  base  drag  as  a  function  of  Mach  number  is  shown  in  Figure  13  and 
is  reproduced  from  Reference  13.  As  shown  in  Figure  13,  reduction  in  base 
drag  due  to  a  cavity  of  about  20%  was  observed  at  subsonic  and  transonic 
speeds.  Whether  or  not  the  base  pressure  is  affected  by  the  presence  of  a 
cavity  in  supersonic  flow  is  not  clear.  The  present  computed  results  show  a 
reduction  in  drag  of  about  25%  at  M  =  .9.  The  effect  of  the  base  cavity  on 
the  base  flow  has  thus  been  correctly  predicted  qua! i tati vely .  Direct  quanti¬ 
tative  comparison  of  the  magnitude  of  the  drag  reduction  has  not  been  made 
since  such  data  is  not  available.  These  results  are  considered  to  be  prelimi¬ 
nary  and  further  computational  effort  is  required. 


VI.  CONCLUDING  REMARKS 

A  computational  predictive  capability  has  been  developed  to  compute  the 
base  flow  for  shell  with  a  base  cavity.  Numerical  computations  have  been  made 
for  the  PXR  projectile  both  with  and  without  the  cavity  at  M  =  .9  and  a  =  0. 

Flowfield  results,  both  with  and  without  the  cavity  have  been  presented. 
The  presence  of  the  cavity  has  been  shown  to  have  a  large  influence  on  the 
near  wake  flow  field.  A  small  recirculation  bubble  is  formed  inside  the 
cavity  which  displaces  the  shear  layer  and  forces  the  rear  stagnation  point 
further  downstream  of  the  base. 

Quantitatively,  the  cavity  has  practically  no  effect  on  the  pressure  drag 
and  the  skin  friction  drag  but  strongly  affects  the  base  drag.  A  reduction  in 
base  drag  and  thus  the  total  drag  of  about  25%  has  been  found  due  to  the  base 
cavity.  The  magnitude  of  this  reduction  may  be  overpredicted  due  to  grid 
resolution  and/or  the  turbulence  model.  Additionally,  oscillations  in  the 
computed  solutions  were  observed  which  are  believed  to  be  numerical  rather 
than  physical.  These  oscillations  seem  to  die  out  when  smaller  and  smaller 
time  steps  are  taken  and  solutions  are  marched  in  time  to  obtain  the  desired 
steady  state  results.  Some  questions  as  to  whether  the  flow  field  in  the  base 
region  is  unsteady  or  not  still  remain.  Future  computations  must  address 
these  issues. 


13. 


S.  N.  B.  Murthy  (Ed.),  "Progress  in  Astronautics  and  Aeronautics:  Aero¬ 
dynamics  of  Base  Combustion ,"  Vot.  40,  AIAA,  New  York,  1976,  pp.  90-92. 
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Figure  2.  Flow  Field  Segmentation  (with  base  cavity) 
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Figure  3.  Model  Geometry 


Figure  4.  Base  Cavity  Configurations  of  Interest 
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Figure  5.  Full  Computational  Grid 
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Figure  6.  Expanded  View  of  the  Grid  Near  the  Model 
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Figure  7.  Expanded  Grid  in  the  Base  Region  (without  Base  Cavity) 
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Figure  8.  Expanded  Grid  in  the  Base  Region  (with  Base  Cavity) 


21 


to 


22 


23 


lit  =  .01 
A  l  =  .005 


60 


s 


2 


<o 

*Q 

m 

in 

»n 


*Q 

K) 

to 

2 


CM 


h 


ir\ 

r^. 

o 


o 

LT\ 

O 


JO 

C\J 

CD 


JO 

CNJ 

o 


I 


CD 
l n 
o 


I 


LTs 

h- 

CD 


I 


24 


Figure  11.  Base  Drag  vs.  Non-dimensional  Time 
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Figure  12.  Total  Drag  vs.  Non-dimensional  Time 


Figure  13.  Base  Pressure  Coefficient  vs.  Mach  Number  under  Different 

Base  Configurations  (Reference  13) 
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